Renormalization and universality of NN interactions in Chiral Quark and 

Soliton Models * 

E. Ruiz Arriola and A. Calle Cordon 

Departamento de Fi'sica Atomica, Molecular y Nuclear, 
Universidad de Granada, E-18071 Granada, Spain 

o 
o 

£NJ ' We use renormalization as a tool to extract universal features of the NN interaction in quark 

and soliton nucleon models, having the same long distance behaviour but different short distance 
components. While fine tuning conditions in the models make difficult to fit NN data, the intro- 
duction of suitable renormalization conditions supresses the short distance sensitivity. Departures 
^ ■ from universality are equivalent to extracting information on the model nucleon structure. 

1. Introduction 

D ' The meson exchange picture has played a key role in the development of Nuclear Physics [1,2] . 

However, the traditional difficulty has been a practical need to rely on short distance information 
which is hardly accessible directly but becomes relevant when nucleons are placed off-shell. From 
a theoretical point of view this is unsatisfactory since one must face uncertainties not necessarily 
■ linked to our deficient knowledge at long distances and which are difficult to quantify. On the 

other hand, the purely field theoretical derivation yields potentials which present short distance 
singularities, thereby generating ambiguities even in the case of the widely used One Boson 
Exchange (OBE) potential. Consider, for instance, the venerable One Pion Exchange (OPE) 
NN — ► NN potential which for r ^ reads 

Vnn jvatM = 7i • t 2 g x ■ a 2 W^{r) + n • r 2 S 12 W^{r) , (1) 

, where the tensor operator S12 — 3<7i • x&2 • x — a\ • a% has been introduced and 



- 

X 



^M = ^%y„(m,r) , W p (r) = ^ &^Y 2 (m w r) . (2) 

Here ^(a;) = e~ x /x and Y 2 (x) — e~ x /x{\ + 3/x + 3/x 2 ) and J^nn = m^g^NN / (2Mjv); 
/^ JVAr /(47r) = 0.07388 for g^NN — 13.08. As we see, the OPE potential presents a 1/r 3 sin- 
gularity, but it can be handled unambiguously mathematically and with successful deuteron 
phenomenology [3]. Nonetheless, the standard way out to avoid the singularities in this and the 
more general OBE case is to implement vertex functions for the meson-baryon-baryon coupling 
(mAB) in the OBE potentials. This correspondins to a folding in coordinate space which in 
momentum space becomes the multiplicative replacement 

V m AB(q) ~> V mA B{q) \jraAB (<? 2 )] 2 • (3) 

where q 2 = q$ — <f is the 4-momentum. Standard choices are to take form factors of the mono- 
pole [1] and exponential [2] parameterizations 



A 

^m°NN(q ) = ^2 _ q2 ' ^rnNNil ) = eX P 



q 2 



A 2 



(4) 



fulfilling the normalization condition T m NN(m 2 ) — 1. Due to an extreme fine-tuning of the 
interaction, mainly in the channel, OBE potential models have traditionally needed a too large 
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guiNN to overcome the mid range attraction implying one of the largest (<~ 40%) SU (3) violations 
known to date. In our recent works [4-9] we discuss how this problem may be circumvented 
with the help of renormalization ideas which upon imposing short distance insensitivity sidestep 
the fine tuning problem and allow natural SU{3) values to be adopted in such a way that form 
factors and heavy mesons play a more marginal role. Contrarily to what one might naively think, 
renormalization reduces the short distance dependence provided, of course, removing the cut-off 
and the imposed renormalization conditions are mutually compatible operations. 

Of course, the extended character of the nucleon as a composite and bound state of three 
quarks has motivated the use of microscopic models of the nucleon to provide an understanding 
of the short range interaction besides describing hadronic spectroscopy; quark or soliton models 
endow the nucleon with its finite size and incorporate basic requirements from the Pauli principle 
at the quark level or as dictated by the equivalent topology [10-13]. While much effort has been 
invested into determining the short range interactions, there is a plethora of models and related 
approximations; it is not obvious what features of the model are being actually tested. In fact, 
NN studies set the most stringent nucleon size oscillator constant value bN = 0.518fm [13] from 
S-waves and deuteron properties which otherwise could be in a wider range bN — 0.4 — 0.6fm. 
This shows that quark models also suffer from a fine tuning problem. In this contribution we 
wish to focus on the common and universal patterns of the various approaches and to show how 
these fine tunings can be reduced to a set of renormalization conditions. 



2. The relevant scales 

From a fundamental point of view the NN interaction should be obtained as a natural solution 
of the 6-q system. However, in order to describe the NN interaction it is far more convenient to 
study two 3-q clusters with nucleon quantum numbers, a procedure also applied in recent lattice 
QCD investigations of the nuclear force [14,15]. NN scattering in the elastic region corresponds 
to resolve distances about the minimal de Broglie wavelength associated to the first inelastic 
pion production threshold, NN — > NNir, and corresponds to take 2Eqm = 2Mn + yielding 
Pcm = \fm^Mjq — 360MeV which means A m ; n ~ 1/vm^M/v = 0.5fm. This scale is smaller than 
In and 27T exchange (TPE) with Compton wavelengths 1.4 and 0.7fm respectively. Other length 
scales in the problem are comparable and even shorter namely 1) Nucleon size, 2) Correlated 
meson exchanges and 3) Quark exchange effects. All these effects arc of similar range and, to 
some extent, redundant. In a quark model the constituent quark mass is related to the Nucleon 
and vector meson masses through M q = Mn/N c — My/2 which for N c = 3 colours gives the 
estimate M q — 310 — 375McV. Exchange effects due to e.g. One-Gluon-Exchange are ~ e - 2M i r 
since they correspond to the probability of finding a quark in the opposite baryon. This follows 
from complete Vector Meson Dominance (for a review see e.g. [16]), which for the isoscalar baryon 
density, psir), and assuming independent particle motion yields 

d 3 xe^(N\p B (x)\N) = 4tt j( drr 2 \ct>(r)\ 2 j (qr) ~ jjrf^ (5) 

suggesting a spectroscopic factor <fi(r) <~ e -Mvr / 2 My / V Airr at large distances. As we have said 
and we will discuss below these effects are somewhat marginal but if they ought to become 
visible they should reflect the correct asymptotic behaviour. In the constituent quark model 
the CM motion can be easily extracted assuming harmonic oscillator wave functions, (f>(r) ~ 
e -6 r /2 [10,11,13] which yield Gaussian form factors falling off much faster than the experimental 
ones. Skyrme models without vector mesons yield instead topological Baryon densities ps(r) ~ 
e -3m»r^7 [12] corresponding to the outer pion cloud contributions which are longest range but 
pressumably yield only a fraction of the radius. In any case quark-exchange looks very much like 
direct vector meson exchange potential which is ~ e~ MvT . 



3. Chiral quark soliton model 

Most high precision NN potentials providing x 2 /DOF <- \ neec [ to incorporate universally 
the One-Pion-Exchange (OPE) potential (including charge symmetry breaking effects) while the 
shorter range is described by many and not so similarly looking interactions [17]. This is probably 
a confirmation that chiral symmetry is spontaneously broken at longer distances than confinement, 
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since hadronization has already taken place. It also suggests that in a quark model aiming at 
describing NN interactions the pion must be effectively included. Chiral quark models accomplish 
this explicitly under the assumption that confinement is not crucial for the binding of tt, N and 
A. Pure quark models including confinement or not have to face in addition the problem of 
recovering the pion from quark-gluon dynamics. In between, hybrid models have become practical 
and popular [10, 11, 13]. As mentioned, all these scales around the confinement scale are mixed 
up. Because these effects are least understood and trigger side effects such as spurious colour Van 
der Waals forces arising from Hidden color singlet states [88]^ states [18, 19] in the (presumably 
doubtful) adiabatic approximation, we will cavalierly ignore the difficulties by remaining in a 
regime where confinement is not expected to play a role and stay with standard chiral quark 
models. 

While both the constituent chiral quark model and the Skyrme soliton model look very dis- 
parate the Chiral Quark Soliton Model embeds both models in the small and the large soliton 
limit respectively El- We analyze the intuitive non-relativistic chiral quark model (NRCQM) ex- 
plicitly and comment on the soliton case where similar patterns emerge. The comparison stresses 
common aspects of the quark soliton model pictures which could be true features of QCD. While 
the long distance universality between both NRCQM and Skyrme soliton model NN calculations 
may appear somewhat surprising this is actually so because in a large N c framework both models 
are just different realizations of the contracted spin-flavour symmetry [23]. 



4. The non-relativistic chiral quark model 

To fix ideas it is instructive to consider the chiral-quark model which corresponds to the Gell- 
Mann-Levy sigma model Lagrangean at the quark level [24] (the non-linear version suggested in 
Ref. [25] will be discussed below), 

C = qffl-g* m ((T + iT 6 T-w))q+± [(d"a) 2 + (d^f] -U(a,n), (6) 

where U(a, tt) = A 2 (er 2 + 7r 2 — v 2 ) 2 /8~ f n m^a is the standard Mexican hat potential implementing 
both spontaneous breaking of chiral symmetry as well as PCAC yielding the Goldberger-Treiman 
relation M q = g^qfir = gaqqfir at the constituent quark level. When this model is interpreted 
from a gradient expansion of the NJL model quarks are regarded as valence quarks whereas 
kinetic meson terms arise from the polarization of the Dirac sea and m 2 = 4M 2 + m 2 , which for 
M g = M^r/3 = My 1 2 yields m a — 650 — 770MeV. In the heavy constituent quarks limit the 
model implies Itt and la exchange potentials, 



2 



'" fp_ 'P)K' -p) 

27r)3 e 
, 1 

;,;r ''"' / (2tt) 3 ~ p 2 +m 2 4tt 



V l7 l(r) = qq r ■ t' I —J—piP-r 

4Af 2 q q J (2tt) 3 p 2 + m 

r J -';(F) : a 2 I e iff,f - ^_^311 (7) 
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m a r 



whence baryon properties can be obtained by solving the Hamiltonian 



N c r 

H = H 

i=l 



P 



2M q 



Mq 



P 2 



^ V(x i -x j ) = — + N c M q + H int , (8) 

i<j 



where the total momentum P — Pi/N c and the intrinsic Hamiltonian have been introduced. 
Due to Galilean invariance the wave function of a moving baryon can be factorized 

* B ( Xl ,...,x N J = . . . , ^ 1 )e lP - R , (9) 



Within the large N c framework the difference corresponds to a saddle point approximation around a trivial 
or non-trivial background. The question which regime is the appropriate one is a dynamical issue [20, 
21]. Likewise, when the soliton is large, quarks are deeply bound and the topological soliton picture of 
Skyrme sets in, giving the appearance of a confined state (where colour Van der Waals forces cannot take 
place) . The soliton of the Spectral Quark model does not allow this interpretation as baryon charge is never 
topological [22]. 



4 
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with R = X)i=i Xi/N c the CM of the cluster and £j = Xi — R/N c intrinsic coordinates, Yli = 0- 
We will assume that this complicated problem has been solved already Ref. [26]. For large N c 

the Hartree mean field approximation b(%i, ■ ■ ■ , %n c ) = TliJii4'oti( x i)Xc might be used [27]). 
For separated hadrons the interaction between quark clusters A and B can be written as sum of 
pairwise interactions which, for elementary irqq and aqq vertices, reads 

Vtofa, ...,x Nc ;y 1 ,..., y Nc ) = £ V^(x, - %) = J -0^ £ v ^ {q yH^y 3 ) . (10) 



Switching to intrinsic coordinates variables xt = ^i+R/2 and yfj = ffj —R/2 with = m — 

where R is the distance between the CM of each cluster, we have 

= %t I ^ iqR ^f L rG'T(q)G^( q r, (11) 

J (27r) d q z + m 2 

2 f Jq-R _ ' 

(2tt) 3 g 2 



^.(i?) = gl qq / ^e^-^-— 5-^(^(9)*, (12 



where the spin-isospin density and scalar densities are given by (e.g. cluster A) 

N c N c 

g a(i) = aE ^ 6 *'" 9 ' PA(?) = ^E ei£i ' 9 > ( 13 ) 

i— 1 i— 1 

respectively. Note that the scalar and Baryon densities as well as the pseudoscalar and axial 
densities coincide unlike the relativistic case. That means that within the approximations one 
should have Ms = My- Thus, the total Hamiltonian is written as 

H = H AM +H BM +V int (R) + 1 ^ r + %-- (14) 

zMt £p 

Galilean invariance implies that inertial masses are Mt — 2N c M q and p — N c M q /2. Intro- 
ducing the two independent cluster complete states HA,int</ ) A,n = MA,n<t>A,n an( i Hb. 
Mb m^B m the two-clusters CM frame unperturbed wave function is just a product 

*^ iBm (1, 2, 3; 4, 5, 6) = <^, n (l, 2, 3; iZ/2)&, (fn (4, 5, 6; -R/2)e iQR , (15) 

where Q is the relative momentum between the two clusters. The above problem is usually 
handled by Resonating Group Methods [10,11,13,28]. We analyze this coupled channel scat- 
tering problem perturbatively where the transition potentials, defined as V An B m - Ak B t {R) = 
(4>A,n4'B,m\Vi n t\4>A,k4>B,i) , have a familiar folding structure which in the case of the pion reads 

VlZ Bm ., AkBl (R) = %§J (2 t)3 ^^/ ^{AnP^Au^B^G^-q)^) . (16) 



5. Long distance limit and the need for renormalization 

At long distances the leading singularities q — im^ and q = im a dominate [29,30]. Using 
that [(TVlp^liV)! 2 is an even function of q we get the structure for the AW — > NN potentials 

VM = 9t n N? [ ^ e ^ \m^pl + C S^(R) + C 2 (-V 2 + ml^HR) + . . . 
qq J (27r) J q- + mi 

= - < —^- + distributions (17) 

An r y ' 

and Eq. |T]) for the OPE contribution. Here, the couplings are given by g a NN — N c g aqq \p(im cr )\ 
and QttN n = N c g A g„ qq \p(im„)\ where g A = (N c + 2)/3 [31]. Assuming \p(im^)\ ~ \p(0)\ = 1 
one has the Goldberger-Treiman relation g A M^ — g^NNf-n at the nucleon level. Thus, at long 
distances finite size effects are represented as an infinite sum of delta functions and derivatives 
thereof. However, any finite truncation will produce a negligible contribution at any non- vanishing 
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distance. In a sense, this result is reminiscent of the Gauss theorem for charged objects with a 
sharp non-overlapping boundary; the interaction is mainly due to the total charge and regardless 
on the density profiles of the system. Only an infinite number of terms may yield a finite size 
effect. Note that the coefficients of the contact interactions are fixed numbers having a meaning 
perturbatively. However, if one tries to play with them to characterize finite resolution effects 
(nucleon size and potential range) in a model independent way non-perturbatively (solving e.g. 
the Schrodinger equation) important restrictions arise. Unlike the S's, the OPE short distance 
1/r 3 singularity is not located in a compact region, i.e. is not killed by taking a finite support 
test function, and contributes to all arbitrarily small distances. Thus, one can effectively drop the 
derivatives of distributions. This simple-minded argument was advanced in Rcf. [32] and explicitly 
verified in momentum space by taking Co and Ci as real counterterms in the Lippmann-Schwinger 
equation in Ref . [29] ; either C2 becomes irrelevant or the scattering amplitude does not converge. 
Therefore, we represent Co as an energy independent boundary condition. The renormalization 
procedure in coordinate space generally corresponds to I) fix some low energy constants such 
as e.g. the scattering length for s- waves, ctQ, at zero energy as an independent variable of the 
potential, 2) integrate in down to an arbitrarily small cut-off radius r c , 3) construct an orthogonal 
finite energy state by matching log-derivatives at r c and 4) integrating out generating a phase-shift 
5o(p) with a prescribed scattering length ao- This prescription is the renormalization condition 
and the procedure of integrating in and out corresponds to evolving along the renormalization 
trajectory. The crucial aspect is that short distance insensitivity is implemented. The tt + a model 
and OBE extensions are analyzed in detail in Refs. [4,5,9] where form factors after renormalization 
are found to be marginal. 



6. Renormalization of Spin-flavour Van der Waals forces 

The non- linear chiral quark model [25] corresponds to take m a — > 00, reducing to just OPE. 
The results for the phase shifts in the lowest partial waves are presented in Fig. [TJ Note the bad 
1 Sq phase. To improve on this the long distance OPE transition potential is taken 

Vab,cd{R) = {tab ■ ? C d){o-ab ■ o-cd[W^]ab ; cd(R) + [Si2]ab,cd[W^] AB ;Cd{R)} , (18) 

where the tensor term is defined as S12 = 3(o~ab • R)(o~cd • R) — <?ab • &cd and 

[W&} AB ;Cd(R) = ^ fwAC ^ BD Y 0t2 (m w R) (19) 

Note that also here there is a 1/r 3 singularity. In this particular form the resulting potential is 
model independent [33] 0- In general, this requires solving a coupled channel problem [34,35] but 
if we are interested in the elastic channel with Tcm = < A = Ma — Mn = 293MeV we may 
take into account the effect of the closed channels as sub-threshold effects in perturbation theory. 
We neglect the exponentially ~ e ~ 2M q r suppressed quark exchange contribution. In obvious 
operator-matrix notation and restricting to the two particle ground |0) = \NN) and excited 
\n) = I ATA), |AiV), |AA) in-going and out-going channels and resolvent Gq^(E) — (E — i/o,fc) _1 
with Ho k — P 2 /(2/j,k) + Ek, we get for the T-matrix 

(T) nm = (V) nm + £ (V) nk G ,k ( V)k, m + 0{V 3 ) , (20) 

k 

with Eq = 2M]\f,Ei t 2 = Mn + Ma and E$ = 2Ma the corresponding thresholds. Thus, separating 
the elastic term k — explicitly from the sum we get the effective potential in the elastic scattering 
channel corresponding to higher pion exchanges, wich, when iterated to second order yields the 
elastic scattering amplitude Too- Specifically, defining the momentum space potential V nm {k' — 
k) = (k',n\V\k,m) = J d 3 RV nm (R)e^ k ~ k ">- R we get 

t7 (V u t?\ \r (V iur f d3q V Qn(k' ~ q)Vno(q - k) / -n 
V 00 (k -k,E) = V 00 (k -*) + Ly ( 27 r)3 E- q y2» n -E n + ° {V ] (21) 



The corresponding couplings are f-wAB = \F 7r AB(im 7r )\ where the transition form factors are defined as 
FnA B {q 2 )x A Tasi XB = {A\G la {q)\B). In the 5(7(4) ® SU C (N C ) quark model [31] and in the chiral limit they 
fulfill .Uaa/Unn = 1/5 and f vNA /f wNN = 3[(N C - l)(N c + 5)/2]*/(iV c + 2). The A — > Nty width in the 
Born approximation yields /^jva/(4"") = 0.324. 



G 
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Fig. 1. Renormalized (eigen) phase shifts for the OPE and A-Born-Oppenheimer potentials as a function of the 
CM np momentum p in the spin singlet 1 S'o (one counterterm) and triplet 3 Si — 3 D\ (three counterterms) channels 
compared to averaged Nijmegen potentials [17]. We take /^]vjv/47t = 0.07388 [17] and /h-Na/ fnNN = 6\/2/5. 



which, expectedly, depends on the energy. Evaluating on-shell at E = Eq + p 2 /2/xq, assuming a 
large splitting p <C V AM a = 600MeV and neglecting the kinetic energy piece in the A A cannel, 
we get the perturbative and local optical potential in coordinate space 

which is the Born-Oppenheimer approximation to second order which generates more complicated 
spin-isospin structures than just OPE including a central force, all of them ~ e" 2m * fl and re- 
sembling TPE. Note that only the intermediate NA state contributes. The above result implies 
an attractive and short distance singular potential since V^ N NA (R) ~ g A /(f 2 R 3 ) and hence 
the potential becomes singular Vnn.nn ~ — Sa/(A/^-R 6 )- Actually, Eq. (|2"2"|) was evaluated in 
the Skyrme soliton model within the Heitler-London approximation, i.e. the product ansatz in 
the coupled channel space [36,37] providing the long sought mid range attraction [12]. We 
reproduce the same results in the quark model calculation. The potential found using Feynman 
graph techniques [39] looks very similar with identical short distance singular behaviour identify- 
ing HaI Qa = fnN&l'Zf-KNN- Note that we leave out background irN scattering which correspond 
to triangle and box TPE diagrams at the quark level. The renormalization procedure as well 
as the necessary counterterms in the general coupled channel singular potentials has been ex- 
plained in much detail in Ref. [32,40]. The results for the phase shifts using Eq. (|22|) in the 
lowest partial waves are depicted in Fig. [TJ In any case the description looks extremely similar 
(including deuteron properties) to the renormalization [41] of more sophisticated field theoretical 
potentials [39]. Convergence is achieved already at r c ~ 0.5fm. 

The multiplicative structures of Eq. (|22|) reflect spin-flavour excitations and remind of the 
analogous Van der Waals forces in atomic systems. They hold literally even after inclusion of form 
factors with folded potentials (although A^^vjv, A^jva and A^-aa are not necessarily identical) 



3 Molecular methods used in the Skyrme model [12,36,37] are replaced by evaluating model form factor yielding 
regularized Meson Exchange potentials [38] where the only remnant of the model is in the meson-form factors. 
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which remove the singularity. This is not equivalent to regularize the effective potential as a whole 
through subtractions. We have checked that form factors after renormalization become marginal 
in agreement with the OBE analysis [9]. 



7. Wigner SU(4) as a long distance symmetry 

If the tensor force component of the qq potential, Eq. 0, is neglected one has invariance under 
the spin-isospin SU (4) group with the quarks in the fundamental 4-dimensional representation, 
q = (u t, u {, d T, d [). In the three quark system we have the spin- flavour states A® A® 
4 = 4a © 20s © 20jv/i © 20m 2 . Due to colour antisymmetry only the symmetric state survives 
which spin-isospin, (S,T), decomposition is 20 s = (|, |) © (§, §) — N A yielding N - A 
degeneracy. Since Ma — Mjv is large at nuclear scales, one might still treat the Nucleon quartet 
N = (p j,p l,n 1,n i) as the fundamental rep. of the old Wigner-Hund SU(4) symmetry which 
implies spin independence, in particular that Vi g (r) = Vs g 1 (r) at all distances suggesting that 
phases Si$ {p) — S3g x (p) in contradiction to data (see e.g. Fig.Q}. The amazing finding of Ref. [6] 
was that assuming identical potentials Vis (f) — VStg^r) for r > r c — > one has 

pcotSi So (p) = , pcot5 3Sl (p) = 1 (23) 

ai So C{p) +T>(p) m Sl C{p)+V(p) 

where the functions A(p), B(p), C{p) and T>(p) are identical in both channels, but the experimen- 
tally different scattering lengths ais = — 23.74fm and a3g 1 — 5.42fm yield quite different phase 
shifts with a fairly good agreement. Thus, Wigner symmetry is broken by very short distance 
effects and hence corresponds to a long distance symmetry (a symmetry broken only by coun- 
terterms). Moreover, large N c [23] suggests that Wigner symmetry holds only for even L, a fact 
verified by phase shift sum rules [6]. In Refs. [7,8] we analyze further the relation to the old Ser- 
ber symmetry which follows from vanishing P- waves in S = 1 channels, showing how old nuclear 
symmetries arc unveiled by coarse graining the NN interaction via the Vjowk framework [42] and 
with testable implications for Skyrme forces in mean field calculations [43]. 

The chiral quark model is supposedly an approximate non-perturbative description, but per- 
turbative gluons may be introduced by standard minimal coupling [13], i</> — > iq) + gfi a ■ \ c a /2 with 
\ c a the AT? — 1 Gell-Mann colour matrices. A source of SU (4) breaking is the contact one gluon 
exchange which yields spin-colour chromo-magnetic interactions (SV,- is the tensor operator), 



4 4rf 1 3 I fij Arrnmj 



i<3 



, 2 



^-sdU*} (24) 



breaking the A — N degeneracy. This short distance terms break also the 1 Sq and 3 Si degeneracy 
of the TV TV system providing an understanding of the long distance character of Wigner symmetry. 
Taking the Wigner symmetric zero energy state and perturbing around it, the previous argument 
suggests that \/as Sl — = 0(Ma — Mm) with a computable coefficient. 



8. Conclusions 

Chiral Quark and Soliton models while quite different in appearance provide some universal 
behaviour regarding NN interactions. If the asymptotic potentials coincide, the main differences 
in describing the scattering data are due to a few low energy constants which in some cases are 
subjected to extreme fine tuning of the model parameters. The success of the model at finite 
energy is mainly reduced to reproducing these low energy parameters. 
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